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INTRODUCTION 


The  recent  demands  on  the  automotive  industry  for  weight¬ 
saving  have  focused  attention  on  the  use  of  polymeric  materials 
reinforced  by  short  lengths  of  reinforcing  fibers.  These  short 
fiber  composite  systems  can  be  molded  to  give  rigid,  lightweight 
structural  components .  The  utilization  of  these  materials  in 
load-bearing  applications  necessarily  directs  attention  to  their 
mechanical  properties  and,  in  particular,  to  the  role  that 
processing  plays  in  establishing  the  mechanical  performance  of 
these  multicomponent  materials . 

Unlike  the  traditional  materials  of  construction,  these 
heterogeneous  materials  exhibit  a  wide  range  of  properties  which 
are  dependent  on  the  initial  fiber/resin/filler  composition  as 
well  as  the  particular  internal  microstructure  developed  during 
fabrication.  This  sensitivity  to  processing  is  manifest  in  both 
the  magnitude  and  directional  dependence  of  the  mechanical, 
thermal,  electrical,  and  transport  properties. 

The  wide  range  of  possible  compositions,  fiber  length 
distributions,  and  fiber  orientation  distributions  that  may  be 
utilized  in  and  generated  during  the  manufacture  of  components 
from  these  materials  precludes  a  total  reliance  on  direct 
laboratory  characterization  of  the  anisotropic  mechanical  prop¬ 
erties.  Consequently,  constitutive  relationships  which  connect 
the  composition  and  processing  dependent  microstructure  to 
mechanical  performance  are  an  important  component  of  the 
technology  for  short  fiber  composite  materials. 

Constitutive  relationships  for  unidirectional  continuous 
fiber  laminates  are  reasonably  well  developed.  For  this  special 
material  system,  the  continuity  of  the  fiber  assures  that  the 
strain  field  parallel  to  the  aligned  fibers  is  essentially 
uniform.  However,  the  stress  and  strain  field  transverse  to 
the  fiber  directions  will  vary  within  the  body.  Consequently, 


the  major  difficulty  in  predicting  the  behavior  of  continuous 
fiber  systems  is  associated  with  predicting  the  Transverse 
Yomg's  Moduli  and  Shear  Moduli.  The  situation  for  short  fiber 
composites  is  considerably  more  complex.  Even  if  all  the  fibers 
were  perfectly  aligned/  the  discontinuities  would  cause  fluctua— 
tions  in  the  stress  and  strain  fields  parallel  to  the  fiber  axes . 

Recently,  Wu  and  McCullough  (Development  in  Composite 
Materials,  Applied  Science  Publishers,  London,  1977)  developed 
improved  variational  treatments  which  provide  general  bounding 
relationships  for  the  effective  elastic  properties  of  a  wide 
variety  of  heterogeneous  materials,  viz,  polycrystalline  metals, 
crystalline  polymers,  continuous  fiber  reinforced  composites, 
short  fiber  composites,  and  particulate  reinforced  polymers. 

Upon  the  specification  of  certain  parameters,  the  boionding 
relationships  yield  families  of  specific  constitutive  relation¬ 
ships  which  contain  all  reported  models  as  special  cases .  As 
would  be  expected,  such  general  constitutive  relationships  are 
somewhat  complex.  Accordingly,  it  is  useful  to  introduce 
simplifications  that  can  yield  reasonable  engineering  estimates 
while  reducing  the  computational  effort  required  to  obtain 
estimates  for  the  anisotropic  Young's  moduli  and  shear  moduli. 

Considerable  simplification  has  been  achieved  through  the 
use  of  an  "Aggregate  Model."  The  important  features  of  the 
Aggregate  Model  are  described  in  a  subsequent  section.  In 
essence,  the  Aggregate  Model  treats  a  short  fiber  composite  as 
a "grainy  metal"  in  which  the  properties  of  the  individual 
grains  are  averaged  over  an  orientation  distribution  to  yield 
the  effective  properties  of  the  bulk  material.  Even  with  this 
simplification,  the  computational  effort  remains  tedious. 
Consequently,  computer  algorithms  have  been  developed  to  facil¬ 
itate  the  prediction  of  properties  from  a  knowledge  of  the 
composition,  state  of  orientation,  and  the  aspect  ratio  of  the 
fibers  of  a  sheet  molding  material.  Two  computational  tools 
have  been  developed:  (1)  an  interactive  FORTRAN  routine  for 
general  use,  and  (2)  a  restricted  routine  for  use  on  hand-held 
calculators  such  as  the  Texas  Instrument  TI-59 . 
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In  the  following  sections^  the  elements  of  material  modeling 

c^v\ 

are  reviewed  and  the  Aggregate  Model  developed  as  an  introduction 

K, 

to  the  notation  used  in  the  subsequent  docTxmentation  of  the 
computational  routine.  Examples  are  provided  to  illustrate  the 
use  of  the  computational  tools  rj  ^ 


BACKGROUND 


Before  proceeding  with  the  development  of  models  for 
predicting  the  behavior  of  sheet  molding  materials  in  terms  of 
composition,  fiber  geometry  and  fiber  orientation,  it  will  be 
useful  to  briefly  review  the  basic  principles  and  notations 
used  for  the  description  of  the  mechanical  behavior  of  materials. 
For  this  purpose,  attention  will  be  first  directed  to  the 
general  characterization  of  the  load-deformation  response  of 
homogeneous  (single  component)  materials.  These  results  will 
be  used  in  subsequent  sections  to  develop  the  effective  load- 
deformation  characteristics  of  heterogeneous  (multicomponent) 
materials. 

Generalized  Materials  Descriptors 

The  load-deformation  response  characteristics  of  an  iso¬ 
tropic  material  (e.g.,  an  amorphous  polymer)  are  traditionally 
described  by  the  Engineering  constants:  The  Young's  modulus,  E; 
the  Shear  modulus,  G;  and  Poisson's  ratio,  v.  The  Young's 
modulus  is  used  to  indicate  the  ability  of  a  material  to  trans¬ 
fer  a  pure  extension  strain  (e)  into  a  pure  tensile  stress  (a) ; 
the  Poisson's  ratio  is  used  to  describe  the  extent  to  which  the 
lateral  dimensions  of  a  body  decrease  in  response  to  a  pure 
extensional  strain.  The  Shear  modulus  is  used  to  describe  the 
ability  of  a  material  to  transfer  a  pure  shear  strain  (y)  into 
a  pure  shear  stress  (t) .  For  isotropic  materials,  these 
descriptors  are  not  independent.  It  can  be  shown  that  if  the 
material  properties  are  the  same  in  all  directions,  then  the 
Engineering  constants  are  related  through  G  =  E/[2(l+v)]. 

Many  material  systems  (e.g.,  drawn  polymers,  continuous 
fiber  reinforced  composites)  exhibit  properties  that  vary  with 
the  direction  in  which  the  load  (or  deformation)  is  applied. 

For  the  current  considerations,  materials  with  orthotropic 
symmetry  are  the  class  of  materials  with  the  lowest  symmetry 
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that  need  be  considered.  Orthotropic  symmetry  is  characteristic 
of  materials  whose  properties  are  equivalent  across  three 
mutually  perpendicular  planes.  For  materials  of  this  symmetry 
class,  it  is  necessary  to  characterize  the  load-deformation 
response  characteristics  along  three  directions  of  the  material 
(e.g.,  the  longitudinal,  transverse,  and  perpendicular  axis). 

The  notation  and  load-deformation  descriptions  for  the  three 
distinct  Young's  moduli  and  the  three  distinct  Shear  moduli 
are  schematically  defined  in  Figure  1. 

In  the  Theory  of  Elasticity  (which  provides  the  theoretical 
basis  for  the  analysis  of  the  mechanical  behavior  of  materials) , 
alternate  sets  of  material  descriptors  are  used:  the  compliance 
array,  S,  and  the  elastic  constant  array,  C.  The  compliance 
array  is  used  to  describe  the  various  deformations  that  result 
from  the  application  of  combined  (or  individual)  loads.  The 
elastic  constant  array  is  used  to  describe  the  various  stresses 
that  result  from  a  general  deformation.  In  the  generalized 
notation,  the  tensile  strain  in  the  "1"  direction  that  results 
from  tensile  stresses  in  the  "1",  "2",  and/or  "3"  direction  are 
given  by 

El  =  +  S^2  ^2  ^13  *^3 


Similarly, 


^2  = 

®21  ^1 

^22 

'^2 

®23  ^^3 

"3  = 

®31  *^1 

^32 

^2 

®33  '^3 

with 

Sji. 

The 

shear 

deformations 

(y) 

are  related  to  the  shear  stresses 

(t)  by  the  relationships 
■^23  ""  ^44  "^23 
^13  ""  ®55  '^13 
^12  ""  ^66  "^12 


FIGURE  1 


Schematic  Definition  of  the  Distinct 
Young's  Moduli  and  Shear  Moduli  for 
Materials  of  Orthotropic  Symmetry 


These  six  relationships  can  be  written  in  compact  matrix  form 
as 


^1 

^11 

^12 

^13 

0 

0 

0 

^1 

"2 

^21 

^22 

^23 

0 

0 

0 

^2 

"3 

^31 

®32 

^33 

0 

0 

0 

X 

^3 

^23 

0 

0 

0 

s,  . 

44 

0 

0 

■^23 

^13 

0 

0 

0 

0 

^55 

0 

1 

■^13 

^12  1 

0 

0 

0 

0 

0 

^66 

^12 

or,  in  symbolic  matrix  notation 

e  =  Sg 

where  e  stands  for  the  1x6  column  vector  of  e's  and  y's;  S  stands 
for  the  6x6  array  of  the  S^^'s;  and  g  stands  for  the  1x6  column 
vector  of  a's  and  t's. 

The  elements  of  the  S  array  are  simply  related  to  the 
Engineering  constants  by  the  relationships  summarized  at  the 
top  of  Table  I . 

In  the  matrix  format,  the  stress,  g,  that  would  result  from 
a  general  deformation,  e,  is  given  by 

g  =  Ce 

Consequently,  the  elastic  constant  array  is  the  "inverse"  Cin  a 
matrix  sense)  of  S,  viz. 


The  relationship  of  the  elements  of  the  C  array  to  the 

Engineering  constants  are  summarized  in  Table  I. 

Most  theoretical  treatments  are  formulated  in  terms  of  the 
elastic  constants,  C. 


TABLE  1 

Relationship  between  Elastic  Constants  and  Engineering  Constants 


Orthotropic  materials 


11  =*  ^ 

S,J  =  -v,2£,'‘ 

£13  =  “Vl3^l’^ 

£44 

=  CjV 

12  =  "V|2£|  * 

£23  =  "*'^23^2  * 

£55 

=  G,-,' 

13  ==  -V,3£i“‘ 

Sij  »  —  VjjSj"' 

£33  -£3"^ 

£46 

=  G.y 

Vy  «  Vj,(£y£,) 

Cji  =  £i£l  ~  (£3/£2)v23]X)  C22  —  ~  ^23 

C|2  *  C21  =*  [£2^12  +  ^3Vi3V23]/>  C23  =*  C32  =  (£3/£|)[£i V23  +  ^55  =  G,3 

Ci3  s=  C3J  =a  £3[Vi2V23  ^  C33  =  £3[l  —  ^^66  =  ^12 

i>’»  =  1  -  2(£3/£JVi2V23V,3  -  Vf3(£3/£|)  -  vi3(£3/£2)  -  V?2(£2/£,) 


Transversely  isotropic  materials 


“I”  axis  unique 

**3**  axis  unique 

£jL,  =*  £i,  £t  “  £2  ^  £3 

£l  «  £3,  £t  ”  £,  =  £2 

Va  *  V,2  «  Vi3.  Vt  «  V23 

Va  =*  Vi3  =  V23,  Vy  =  V12 

”  ^12  =“  G|3,  Gy  G23 

£t/[2(I  +  Vy)]  Gy^  =s  Gi3  =»  G23,  Gy  =s  Gj2 

Isotropic  materials 

-  £x/[2(l  +  yj)] 


E  ^  El  ^  Ez  ^  Ei 

V  *  V,2  =  Vji  =  V,3  =  V31  =r  V23  =  V32 

G  ^  Gi2  ^  G,3  =»  G23  =  £/[2(I  +  v)] 


(From:  Anisotropic  Elastic  Behavior  of  Crystalline  Polymers, 

R.  L.  McCullough,  Treatise  on  Materials  Science  and  Technology, 
lOB,  p.  453-540,  Academic  Press,  New  York,  1977) 


Additional  symmetry  of  the  material  results  in  certain  special 
conditions  on  the  elastic  constants,  C  (as  well  as  the  compliance 
constants,  S,  and  the  engineering  constants) .  Transversely, 
isotropic  symmetry  is  of  particular  importance,  fibers  (e.g., 
graphite,  Kevlar)  frequently  exhibit  transversely  isotropic 
response.  Under  conditions  of  transverse  isotropy,  the  structure 
of  the  orthotropic  array  is  preserved;  however,  certain  special 
relationships  between  the  material  descriptors  are  imposed. 
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These  relationships  are  summarized  in  Figure  2  and  in  Table  I. 

Isotropic  symmetry  is  exhibited  by  many  (undrawn)  polymer 

systems.  Again,  the  structure  of  the  orthotropic  array  is 

preserved  with  the  special  conditions  =  C22  “  ^33' 

C,.  =  =  C^r  =  Consequently,  no  unique  material 

44  55  66  11  12 

axes  exist  for  isotropic  materials  (i.e.,  the  properties  are 
the  same  in  all  directions) .  The  results  of  the  symmetry 
conditions  for  isotropic  materials  are  summarized  in  Figure  2 
and  Table  I . 

These  descriptions  may  be  applied  to  describe  the  mechanical 
properties  of  the  individual  components  of  a  sheet  molding 
compound  as  well  as  the  overall  properties  of  the  sheet  molding 
material. 

Orientation  Dependence 

In  the  preceding  discussion  of  the  characterization  of 
material  descriptors,  the  simplifications  resulting  from 
material  symmetry  were  emphasized  by  requiring  that  the  unique 
axis  associated  with  the  symmetry  class  of  the  material  were 
coincident  with  the  loading  (or  deformation)  directions  imposed 
on  the  body  of  the  material.  This  situation  is  rarely  encoun¬ 
tered.  For  example,  fiber-reinforced  laminates  (e.g.,  tire 
plies)  are  frequently  oriented  at  an  angle  of  45°  with  respect 
to  the  fiber  direction.  Fortunately,  once  the  material  system 
has  been  characterized  along  the  unique  axes,  the  material 
response  in  any  arbitrary  direction  can  be  accurately  predicted 
through  the  relationships  summarized  in  Figure  3. 

These  relationships  can  be  used  to  characterize  the  load- 
deformation  response  characteristics  of  a  material  A((}))  at  any 
arbitrary  direction,  ({> ,  in  terms  of  the  associated  descriptors 
associated  with  the  symmetry  axes  of  the  material,  A. 
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FIGURE  2 


Components  of  the  Elastic  Constant  Array 
for  Materials  with  Orthotropic,  Transversely 
Isotropic,  and  Isotropic  Materials  Symmetry 
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ORTHOTROPIC  MATERIALS 
3  2 


I 


TRANSVERSLY 

ISOTROPIC 

MATERIALS 


k, 


SPECIAL  CONDITIONS 

^22*^33  ^11  "^22 

C|2  =  C|3  C23=C|3 

C44=  Y  ^^22*^23^  ^44’  ^55 
^55=^66  Ce6= 


ISOTROPIC 

MATERIALS 

3 


SPECIAL  CONDITIONS 

C||  =  C22  =  C33 
C12  '  ^13  ■  ^23  j 

C44=  C55  =  Ce<3=  ^(C||-C|2) 
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3,3' 


A||  «  *  BqCOSA^ 

A|2  -  Azi  •  Bj»-  B,cos44 
^le  *  ^«r  *  '^B7sin2^)4>  BgSin4^ 
A22  *  B1-B7COS24  ■¥  B9C0344 
A26  ■  Afij*  ■^{B7sin24)*B3im4^ 
^66  *  Bg-B0CO*4^ 


(W|l  •  CO94  <**|2  “  *'^13  * 

4(21  •**s*n^  **ZZ’^  cos^  ^^23*  ®  I 
•»3I*0  <1132*0  <i»33”  1/ 

3}  *  '^(3A|| 3A22'*’2A|2'*‘4A05} 
B7  »  A||-A22) 

Be  *  i«  A,,  ♦  A22-2Ai2.4A6e) 
82  *  i  (  All  ♦  A22 ‘*'6A|2<>4Ae0) 
Be  •  8  *  A||  ♦  A22*2A|2-*’4A66) 


Reduction  in  the  transformation  relationships  for  an  orthotropic  material  in  a  state 
of  plane  stress  or  strain.  The  rotation  transformation  between  the  material  axes  (ej ,  e2,  ^3)  and 
arbitrary  load  (or  deformation)  axes  (I,',  ^2',  e^)  reduce  to  a  simple  rotation  around  the  common 
^3  axis.  Replacement  of  the  Aij  by  elastic  constants  C,j  yields  expressions  for  the  transformed 
elastic  constants  «  Cy);  replacement  of  the  /4y  by  the  compliance  constants  Sy,  yields 
expressions  for  the  transformed  compliance  constants  {Aij  =:  Sfjlmimj)\  m*  =  I  for  fc  =  2, 

or  3;  =  2  for  k  ^  4,5,  or  6. 


FIGURE  3 

(From:  Anisotropic  Elastic  Behavior  of  Crystalline  Polymers, 

R.  L.  McCullough,  ob.  citO 


MATERIALS  MODELING 


In  this  section,  procedures  will  be  developed  for  predicting 
the  behavior  of  sheet  molding  materials  comprised  of  short 
reinforcing. fibers,  particulate  fillers,  and  a  polymer  phase. 

The  approach  proceeds  by  developing  an  aggregate  model  to  account 
for  orientation  effects.  The  basic  element  of  the  aggregate 
model  is  taken  to  be  an  arbitrary  "micro-region."  The  properties 
of  the  micro-region  are  subsequently  predicted  by  associating 
the  micro-region  with  a  micro-laminate  of  aligned  fibers .  The 
combination  of  these  results  into  a  computational  format  for 
predicting  the  properties  of  sheet  molding  materials  is  summar¬ 
ized  in  the  final  section. 

Aggregate  Model 

Under  the  Aggregate  Model,  a  sheet  molding  material  is 
viewed  as  partitioned  into  micro-regions  as  illustrated  in 
Figure  4.  As  in  the  case  of  "grainy  metals,"  each  micro-region 
is  treated  as  an  apparent  homogeneous  (but  anisotropic)  material 
which  may  be  described  as  an  array  of  elastic  constants,  C. 

These  micro-regions  may  assume  a  variety  of  orientations  with 
respect  to  the  external  axes  of  the  macroscopic  body  of  material . 
Accordingly,  the  load-deformation  characteristics  along  the  "1", 
"2",  and  "3"  axes  of  the  bulk  material  are  dependent  upon  the 
relative  orientation  of  the  unique  "1"',  "2'",  "3"'  material 
axes  of  the  micro-region.  Consequently,  the  transformations 
given  in  Figure  3  must  be  applied  to  each  micro-region  to 
obtain  the  appropriate  contribution  of  the  individual  regions 
to  the  overall  load-deformation  response  of  the  bulk  material. 

The  fraction  of  micro-regions  whose  unique  material  axis 
"1'"  makes  an  angle  (j)  with  respect  to  the  external  body  axis  "1" 
may  be  specified  in  terms  of  an  orientation  distribution 
function  n((f))  .  This  function  gives  the  relative  number  of 
micro-regions  whose  unique  axes  are  parallel  and  make  an  angle 


FIGURE  4 


Schematic  Definition  of 
the  Aggregate  Model 
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(f)  with  respect  to  the  external  axis.  Thus,  if  all  the  micro¬ 
regions  were  aligned  along  the  body  axis,  n((()=0)  =  1  and 
n(<t)7^0)  =  0.  Alternately,  if  the  micro-regions  were  uniformly 
distributed  in  all  directions,  n((j))  =  a  constant  for  all  values 
of  <p . 

It  should  be  emphasized  that  the  orientation  function, 
n((}))  ,  does  not  provide  for  an  "out-of-plane  tilting"  of  the 
micro-region.  This  type  of  orientation  function  is  called  a 
"Planar"  distribution.  Planar  distributions  are  appropriate 
to  sheet  molding  materials  that  are  fabricated  under  conditions 
which  maintain  such  planarity.  Materials  formed  by  injection 
molding  may  exhibit  "out-of-plane  tilting"  and  therefore  will 
require  a  different  form  of  an  orientation  distribution 
function. 

The  planar  distribution  function  has  the  following  important 
features : 

n(<J))  =  n(-({i) 

n((j))  =  n(+(t)) 

ir/2 

n((p)d(p  =  1 

0 

The  net  contribution  of  the  various  micro-regions  to  the 
overall  load-deformation  response  is  given  by  averaging  the 
relationships  given  in  Figure  3  over  the  distribution;  viz 

* 

<A>  =  A(4))n((j>)di})  ....1 

It  is  useful  to  introduce  certain  orientation  parameters, 

"f"  and  "g"  defined  as 

2 

f  =  2  <cos  (j)>  -1  ....  2 

g  =  (1/5)  (8  <cos^<j)>  -3)  -  2b 

with  it/2 

<cos'^(})>  =  cos^(})n  (()))  d({) 


0 


These  orientation  parameters  are  constructed  to  provide  a 
convenient  scale  for  characterizing  the  state  of  orientation 
of  a  system.  Thus,  for  f  =  g  =  0,  the  micro-regions  are  randomly 
distributed  within  the  "1-2"  plane  of  the  bulk  material.  For 
f  =  g  =  1,  the  micro-regions  are  perfectly  aligned  along  the 
"1"  axis  of  the  bulk  material.  Values  of  f  and  g  between  0  and 
1  represent  intermediate  states  of  orientation.  These  features 
of  the  planar  orientation  are  summarized  in  Figure  5. 

The  results  obtained  by  averaging  under  a  planar  orientation 
distributions  are  summarized  in  Table  II  in  terms  of  the  orien¬ 
tation  parameters  f  and  g.  As  before,  the  elastic  constant 
array,  C,  and  the  compliance  array,  S,  undergo  the  same  trans¬ 
formations  so  that  the  results  for  orientation  averaging  can  be 
generalized  for  an  arbitrary  material  descriptor,  <A>  and  base 
descriptor,  A.  Thus  for  A  ^  C,  the  relationships  of  Table  II 
yield  the  orientation  average  of  the  elastic  constant  array; 
for  A  -»■  S,  the  orientation  average  of  the  compliance  array  is 
obtained.  The  factor  B  is  introduced  to  account  for  the  con¬ 
traction  of  the  compliance  and  elastic  constants  to  second 
order  tensors.  For  A  C,  B  =  1;  for  A  ^  S,  B  =  4. 

It  can  be  shown  that  for  reasonable  forms  for  the  planar 
orientation  distribution  (e.g.,  n  ((f))  =  kcos^cfi)  that  the 
orientation  parameters  f  and  g  are  related: 

g  =  2f  (7-2f)/[5(4-2f)  ]  _  3 

so  that  only  one  parameter  is  required  to  describe  the  state 
of  orientation.  The  orientation  parameter,  "f",  is  related  to 
the  "root-mean-square"  orientation  angle 

<|)rms  =  cos"^/(l+f)/2  . 

Hence  for  random  distributions,  f  =  0  and  =  45°.  For 

iTins 

perfect  alignment,  f  =  1,  and  <{>„__  =  0°.  For  "slight"  orien- 

X  XILS 

tation,  f  =  0.25,  d)  =  40°;  for  "moderate"  orientation, 

rms 


FIGURE  5 


Summary  of  the  Features  of  a 
Planar  Orientation  Distribution 


TABLE  II 

Averaged  Properties  of  a 
Planar  Orientation  Distribution 
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for 


"significant"  orientation. 


The  aggregate  model  serves  as  a  reasonable  means  of 

introducing  relationships  to  account  for  the  particular  internal 

state  of  orientation  of  the  representative  micro- regions .  The 

next  task  is  to  identify  the  nature  of  the  arbitrary  micro- 

* 

regions  and  relate  the  properties  of  the  micro-region,  C  ,  to 
the  composition  and  fiber  geometry  of  the  sheet  molding  material 


Micro-Laminate  Model 

The  preceding  treatment  of  the  aggregate  model  was  based 
upon  an  arbitrary  specification  of  a  "micro-region."  Indeed, 
the  elements  of  the  elastic  constant  array,  C,  that  is  used  to 
describe  the  mechanical  behavior  of  the  micro-region  could  be 
treated  as  adjustable  parameters  for  curve  fitting.  Since  as 
many  as  nine  independent  parameters  could  be  required,  this 
empirical  approach  does  not  appear  to  be  fruitful .  In  this 
section,  the  characteristic  micro-region  of  the  aggregate  model 
is  treated  as  a  micro-laminate  of  perfectly  aligned  fibers. 
Accordingly,  the  domain  of  a  micro-region  is  specified  as  that 
portion  of  the  material  that  can  be  partitioned  into  volume 
elements  in  which  the  fibers  within  the  region  are  all  aligned 
parallel  to  an  axis  that  makes  an  angle  (j)  with  the  external 
body  axis.  The  identification  of  such  a  region  specifies  the 
effective  aspect  ratio  of  the  fiber.  Thus  for  a  sheet  molding 
material  (e.g.,  SMC-25)  comprised  of  relatively  long  and 
straight  collections  of  fibers,  the  aspect  ratio  would  be  large 
(e.g.,  a  >  250).  Alternately,  for  a  sheet  molding  material 
in  which  the  fiber  buridles  are  "swirled"  (e.g.,  SMC-65)  such 
that  arc  length  is  relatively  short,  the  effective  aspect 
ratio  could  be  as  low  as  a  =  1  ->•  10. 

Under  this  view  of  the  micro-regions ,  the  properties  can 
be  obtained  from  appropriate  models  for  laminates .  The  general 
Wu-McCullough  relationship  can  be  specialized  to  this  purpose. 


The  general  relationship  is  of  the  form 

C*  =  Cq  +  (m"^  +  _  4 

where  C*  is  the  6x6  array  of  elastic  constants  for  the  micro¬ 
laminate.  The  term  is  a  6x6  array  of  elastic  constant  for  a 
"reference"  material.  The  term  E^  is  a  6x6  array  which  takes 

into  account  correlation  effects.  The  elements  of  E  are 

asO 

dependent  upon  the  aspect  ratio,  a,  and  certain  elements  of 

C  .  The  elements  of  E  are  summarized  in  Table  III. 
zo  '  :so 

The  quantity  M  is  a  6x6  array  that  is  dependent  upon  the 
composition  and  the  properties  of  the  components  compensated 
for  correlation  effects,  viz. 


M  =  Zv. m.  ....  5 

»  ISl  ' 


where  is  the  volume  fraction  of  the  i'th  component;  the  term 
m^  is  a  6x6  array  of  the  properties  of  component  "i"  compensated 
for  correlations  through  the  following  relationships 


m. 

a:l 


....  6  a 


li 


=  9i  - 


6b 


The  term  is  the  6x6  array  of  elastic  constants  for  the  i'th 
component 

The  versatility  of  the  Wu-McCullough  relationship  is 

manifest  through  the  term  C  .  Assigning  a  reference  material 

"zero"  rigidity  (C  =  0)  yields  the  classic  Reuss  model; 

assigning  the  reference  material  an  infinite  rigidity  (C^  = 

yields  the  Voigt  model.  If  the  reference  material  is  taken  as 

the  resin  phase  (C  =  C  .  ) ,  the  "best  lower  bounds"  are 

^resin 

obtained.  If  the  fiber  phase  is  selected  as  the  reference 
phase  (C^  =  ,  the  "best  upper  bounds"  are  obtained. 

Usually,  these  bounds  are  too  far  apart  to  provide  useful 
predictions.  For  the  case,  Q*  =  the  "self-consistent" 

field  models  are  obtained. 


TABLE  III 

Elements  of  E 

zo 

n  Z  Z  0  0  0 

Z  k+u  k-u  0  0  0 

Z  k-u  k+u  000 

0  0  0  4u  0  0 

0  0  0  0  4m  0 

0  0  0  0  0  4m 

n  =  4ah^(a)  -43h2(a) 

Z  =  2ah^(a) 
k  =  ah^Ca)  -  3hj^(a) 
u  =  Jsah^Ca)  -  3h^{a) 
m  =  2ah^(a)  -  3[%h^(a)  +  h2(a)] 

“  “  ^^22  ~  ^44^/ ^^*^22^44^ 
s  =  V<4C5^) 

h2(a)  =  1  -  h^(a) 

h^(a)  =  h[l  -  h^Ca)  -  hg(a)] 
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TABLE  III  (con't) 

For  0  _<  a  <  1 

2  2,,,  2, 
y  =  a  /(1-a  ) 

h^(a)  =  y^{[(l/y)  +  y]tan“^(l/y)  -  1} 

h3(a)  =  y^{  t(l+y^)/2y^]  +  1  [hd/y)^  -  (1-y)  -  (3/2)y]  tan"^  (1/y)  ] } 
h5(a)  =  (l+y2)  [(3-y2)/(l+y^)  -  ( 3/2) ytan"^ (1/y) ) ] 

For  a  =  1 

hj^(a)  =  2/3 
h^Ca)  =  8/15 
h3(a)  =  1/5 

For  1  <  a  <  " 

=  (a^-l)/a^ 

Z  =  {£n[  (l+x)/(l-x)  ]  }/x 
hj^(a)  =  [1  -J2(l-x)Z]/x 

h3(a)  =  [(3/2)  -h  +  3«(x^+2x-3)  Z]/x^ 

hgCa)  =  [(l-x^)/x^]^{[(3-2x)/2(l-x)  ]  -  (3/4)Z)} 


It  has  been  shown  that  the  properties  for  a  wide-range  of 

fiber  reinforced  resin  systems  can  be  accurately  predicted  by 

assigning  values  to  the  array  that  correspond  to  a  material 

(of  equivalent  composition  and  concentration)  reinforced  by 

spheres  (aspect  ratio  -  1)  rather  thar  fibers.  It  was  proposed 

that  standard  samples  of  glass  bead  reinforced  resins  be 

prepared  and  characterized  over  a  range  of  volume  fractions  to 

provide  the  data  necessary  for  constructing  the  array. 

These  experimentally  determined  values  for  C  could  be  used  in 

o 

conjunction  with  Eq.  4  to  predict  the  properties  of  fiber 
reinforced  systems. 

Recently,  it  was  shown  that  the  experimental  determination 

of  is  not  required.  A  model  for  particulate  systems  has  been 

developed  which  accurately  predicts  the  behavior  of  a  wide 

variety  of  particulate  filled  systems  over  the  sensible  range 

(v  <  0.8)  of  concentration  of  filler,  v  . 

P  -  '  p 

The  "S-Mixing  Rule"  model  for  particulate  systems  is  of 

the  form 


Ip  " 


V  S_  + 

r»Lo 


where  v  and  v  are  the  respective  volume  fractions  of  resin 
^  P  —  I 

and  filler  particles.  The  quantity  S-  =  C..  where  C,  is 

:;IiO  aLo  aLo 

obtained  from  Eq.  4  with  C  =  C  .  and  the  aspect  ratio  of 
the  E^  term  taken  as  a  =  1.  Similarly,  the  quantity 
is  obtained  from  Eq.  4  with  C  =  C, ...  and  a  =  1.  The 
relationships  given  in  Table  I  are  used  to  obtain  the  Young's 
modulus.  Shear  modulus,  and  Poisson's  ratio  from  the  computed 
values  of  S^  for  the  particulate  system. 

The  demonstrated  success  of  the  "S-Mixing  Rule"  provides  a 
convenient  means  for  generating  the  appropriate  values  for  the 
parameters  of  the  reference  phase,  C^,  as  required  by  Eq.  4. 

For  a  two-component  fiber/resin  system,  the  effective 
mechanical  properties  of  the  material  may  be  predicted  by  the 
following  procedures: 


(i)  For  the  current  volume  fraction  of  fiber  and  resin 
the  appropriate  values  of  are  obtained  by  the 
application  of  the  "S-mixing  rule"  (Eq.  5)  for  a 
system  comprised  of  the  equivalent  volume  fraction 
of  resin  and  particles  (a  =  1)  with  the  particles 
assigned  the  properties  of  the  fiber. 

(ii)  The  values  obtained  for  are  used  in  Eq.  4  along 
with  the  effective  aspect  ratio  of  the  fiber  to 
predict  the  properties  of  the  micro-laminate,  Q* . 

(iii)  The  properties  of  the  micro- laminate  are  subjected 
to  the  orientation  averaging  prescribed  in  Table  II 
for  a  specified  state  of  orientation,  f. 

(iv)  The  averaged  values  of  C  for  the  sheet  molding 
material  are  converted  to  the  compliance  array,  S, 
and  subsequently  to  the  Engineering  constants  via 
the  relationships  given  in  Table  I . 

These  procedures  can  be  extended  to  resin/fiber/filler 
systems  by  introducing  the  notion  of  a  "surrogate"  matrix.  In 
this  approach,  the  resin  and  particulate  filler  system  are 
viewed  as  a  matrix  phase.  The  properties  of  the  matrix  phase 
may  be  predicted  by  the  application  of  the  S-mixing  rule  for  a 
particulate  system  with  the  apparent  volume  fractions  of  v^ 

and  v'  for  the  resin  and  filler.  These  apparent  volume 

p 

fractions  are  related  to  the  true  volume  fractions,  v^  and  v^, 
through  the  relationships 


v'  =  V  /(v  +  V  ) 

r  r  r  p 


V'  =  Vp/(v^  +  Vp) 


so  that  v'  +  v'  =  1. 

The  resulting  properties  for  the  isolated  resin/filler 

systems  are  used  to  represent  the  behavior  of  a  surrogate 

matrix  material  with  properties  as  predicted  by  Eq.  5  with 

concentrations  v'  and  v' .  At  this  point,  the  computation 

r  p 


follows  steps  (i)  through  (iv)  for  a  two-component  fiber/matrix 
system  with  the  computed  values  for  the  surrogate  matrix 
assuming  the  role  of  the  resin  phase.  The  volume  fraction  of 
the  surrogate  matrix  phase  is  given  by  ^  ^p  ~ 

where  v^  is  the  volume  fraction  of  the  fiber  phase. 


FORTRAN  PROGRAM 


31 


Introduction 

The  preceding  sections  have  been  concerned  with  developing 
a  model  for  predicting  the  Engineering  properties  of  sheet 
molding  compounds .  These  ideas  have  been  implemented  in  the 
FORTRAN  program  SMC-3.  Use  of  SMC-3  is  described  in  this 
section.  Two  examples  using  SMC-3  and  a  program  listing  follow. 
Finally,  a  few  cautions  which  need  to  be  noted  when  implementing 
SMC-3  on  a  computing  system  other  than  the  DEC-system  10  for  which 
this  version  was  written  are  discussed. 

Examples  and  use  of  SMC-3  will  be  illustrated  by  specific 
examples .  In  general  the  execution  of  the  program  requires 
the  following  input  data: 

1)  number  of  components 

2)  properties  of  constitutive  phases 

3)  voliime  or  weight  fractions  of  components 

4)  effective  aspect  ratio  of  fibers 

5)  orientation  of  fibers 

To  facilitate  the  use  of  SMC-3,  the  Engineering  constants 
(in  psi)  for  polyester  resin.  E-glass  fibers,  and  calcixim 
carbonate  filler  have  been  stored  internally  in  the  program. 

If  these  properties  are  to  be  used,  no  data  regarding  the 
engineering  properties  of  the  constitutive  phases  need  be 
entered.  If  the  properties  of  one  or  all  the  phases  are  to  be 
changed,  this  can  be  accomplished  by  the  user  during  program 
execution.  For  an  isotropic  phase,  e.g.,  resin,  filler  and 
some  fibers.  Young's  modulus,  the  Shear  modulus  and  Poisson's 
ratio  will  need  to  be  entered.  For  transversely  isotropic 
fibers  (e.g.,  Kevlar)  two  Young's  moduli,  two  Shear  moduli  and 
two  Poisson's  ratios  will  need  to  be  entered.  Example  1 
demonstrates  how  properties  for  isotropic  phases  are  changed. 
Example  2a  shows  how  the  properties  are  changed  for  a  transversely 
isotropic  fiber. 

The  composition  can  be  entered  either  as  volume  fractions 
or  weight  fractions.  If  volume  fractions  are  used,  the  data 
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are  entered  directly.  Input  using  weight  fractions  also  requires 
that  the  density  of  each  phase  must  be  entered.  In  examples  1 
and  2,  composition  is  input  as  weight  and  volume  fractions, 
respectively. 

The  effective  aspect  ratio  is  entered  directly  upon  request. 
The  aspect  ratio  ranges  from  one  to  infinity.  An  aspect  ratio 
of  one  corresponds  to  a  spherical  inclusion  while  an  infinite 
aspect  ratio  corresponds  to  a  continuous  fiber. 

Effect  of  fiber  orientation  on  the  Engineering  constants 
can  be  investigated  in  either  of  two  modes.  First,  the  Engi¬ 
neering  constants  can  be  calculated  for  a  single  user  selected 
fiber  distribution.  Orientation  is  specified  by  Herman's 
orientation  function  f.  For  random  distribution  f  =  0.  For 
perfectly  aligned  fibers  f  =  1. 

Slightly  oriented  systems  can  be  represented  by  "f"  values 
in  the  range  0.2  to  0.3.  Moderately  oriented  systems  can  be 
represented  by  "f"  values  of  ~0.5.  Significantly  oriented 
systems  can  be  represented  by  "f"  values  of  0.6  to  0.8 

Second,  a  range  of  fiber  orientations  can  be  scanned.  In 
this  case  the  user  specifies  the  initial  and  final  values  of 
the  orientation  parameter  as  well  as  the  incremental  step  size. 

For  some  input  data  it  is  possible  to  perform  consistency 
checks.  SMC-3  provides  two  such  checks.  The  first  is  for 
composition  fractions.  If  the  fractions  do  not  sum  to  unity, 
the  user  is  instructed  to  re-enter  the  data. 

The  second  consistency  check  is  performed  when  one  of  the 
standard  phases  is  replaced  by  an  isotropic  phase.  In  this 
case  there  are  only  two  independent  material  descriptors,  and 
it  is  possible  to  determine  whether  the  input  Young's  modulus. 
Shear  modulus  and  Poisson's  ratio  are  self-consistent.  The 
results  of  this  check  are  as  follows : 

1)  If  input  data  are  self-consistent,  program  execution 
continues . 


2)  If  the  input  Poisson's  ratio  is  inconsistent,  an 
internally  determined  value  is  assigned  and  the  user 
is  informed  of  the  change  and  asked  for  confirmation. 

3)  If  the  input  Young's  and  Shear  moduli  require  a 
negative  Poisson's  ratio  or  one  larger  than  0.5,  the 
user  is  instructed  to  enter  new  moduli  and  Poisson's 
ratio. 

Output  from  SMC-3  consists  of  two  parts.  The  first  is  a 
summary  of  the  input  data.  Included  are  the  phase  properties, 
composition  as  volume  fractions,  and  the  effective  aspect  ratio. 
Secondly,  a  tabular  summary  of  the  predicted  material  properties 
at  specified  orientations  is  printed.  The  table  contains  the 
predicted  longitudinal,  transverse,  and  perpendicular  Young's 
moduli,  the  "1,  2",  "1,  3"  and  "2,  3"  Shear  moduli  and  the 
"1,  2",  "1,  3"  and  "2,  3"  Poisson's  ratios. 


EXAMPLES 


EXAMPLE  1 


The  purpose  of  this  example  is  to  illustrate  the  use 
of  SMC3  for  a  three  component  system.  The  input  data  for 
the  filler  will  be  altered.  An  intentional  error  will  be 
introduced  to  illustrate  the  operation  of  the  self-consis- 
tancy  check.  After  illustrating  the  checking  feature,  the 
data  will  be  reset  to  the  standard  values. 

□  Data  inputs  by  the  user  will  be  indicated  by  the  symbol  Q . 

□  Comments  are  given  in  italics. 


ENTER  NUMBER  OF  COMPONENTS? 

FOR  RESIN/FIBER  SYSTEM  NUMB  =  2 

FOR  RESIN/FIBER/FILLER  SYSTEM  NUMB  =  3 

03 

TYPICAL  PROPERTIES  FOR  A  THREE  COMPONENT 
POLYESTER/E-GLASS/CALCIUM  CARBONATE  SYSTEM  ARE? 


YOUNG'S  MODULUS  SHEAR  MODULUS  POISSON'S  RATIO 


RESIN  .5100E+06  .1960E+06  *301 
FIBER  ♦1050E+08  ♦3940E+07  .333 
FILLER  .6930E+07  ♦2620E+07  .323 


IF  THESE  VALUES  ARE  ACCEPTABLE »  ENTER  A  "0*  (ZERO) 

IF  YOU  WISH  TO  USE  SIGNIFICANTLY  DIFFERENT  VALUES  ENTER  A  "1*  (ONE) 

□For  illustrative  purposes^  the  user  elects  to  reduce  the 
Young's  Modulus  of  the  filler  from  6.9Z  X  10^  to  4  X  10^ 
while  maintaining  all  other  values  for  the  properties . 

This  will  give  inconsistent  values  for  an  isotropic 
material . 

Oi 

IF  RESIN  PROPERTIES  ARE  ACCEPTABLE  ENTER  A  *0“  (ZERO) 

IF  THEY  ARE  TO  BE  ALTERED  ENTER  A  "I"  (ONE) 

Oo 


IF  FIBER  PROPERTIES  ARE  ACCEPTABLE  ENTER  A  "0"  (ZERO) 

IF  THEY  ARE  TO  BE  ALTERED  ENTER  A  *1“  (ONE) 

IF  FILLER  PROPERTIES  ARE  ACCEPTABLE  ENTER  A  "0*  (ZERO) 

IF  THEY  ARE  TO  BE  ALTERED  ENTER  A  "I*  (ONE) 

ENTER  YOUNG ^S  MODULUS  FOR  FILLER 

ENTER  SHEAR  MODULUS  FOR  FILLER 
Q2*62E6 

ENTER  POISSON'S  RATIO  FOR  FILLER 

0*323 

*  *  INPUT  ERROR  *  * 

FOR  INPUT  VALUES  OF  THE  YOUNG'S  MODULUS  AND  SHEAR  MODULUS 
POISSON'S  RATIO  WOULD  BE  NEGATIVE  OR  GREATER  THAN  0.5 

RE-ENTER  DATA 

QReturn  to  standard  values. 


ENTER  YOUNG'S  MODULUS  FOR  FILLER 
06*93E6 

ENTER  SHEAR  MODULUS  FOR  FILLER 
0*262E7 

ENTER  POISSON'S  RATIO  FOR  FILLER 
0*323 

THE  CURRENT  SET  OF  PROPERTIES  FOR  THE  THREE 
COMPONENT  SYSTEM  AREJ 


YOUNG'S  MODULUS  SHEAR  MODULUS  POISSON'S  RATIO 


RESIN  .5100E+06 
FIBER  .1050E+08 
FILLER  .6930E+07 


1960Ef06 

♦  301 

3940E+07 

.333 

2620E+07 

.323 

Oo 


IF  THESE  VALUES  ARE  ACCEPTABLE  f  ENTER  A  ■‘O"  (ZERO) 

IF  YOU  WISH  TO  USE  SIGNIFICANTLY  DIFFERENT  VALUES  ENTER  A  'I*  (ONE) 


□  fraation  variables  wilt  be  used.  TJie  use  of  weight 

fraotion  variables  requires  input  for  the  density  (or  speoifio 
gravity)  for  each  component. 


IF  THE  COMPOSITION  IS  TO  BE  ENTERED  AS  OOLUME  FRACTIONS f 
ENTER  A  “0"  (ZERO) 

IF  THE  COMPOSITION  IS  TO  BE  ENTERED  AS  WEIGHT  FRACTIONS » 
ENTER  A  "I*  (ONE) 


WEIGHT  FRACTION  OF  RESIN 

DENSITY  OF  RESIN 

©.t«2 

WEIGHT  FRACTION  OF  FIBER 

©♦25 

DENSITY  OF  FIBER 
©2*55 


WEIGHT  FRACTION  OF  FILLER 
©♦418 

DENSITY  OF  FILLER 
©2*40 

^  ENTER  ASPECT  RATIO  OF  THE  FIBER 
©500^ 

\3This  value  for  an  aspect  ratio  is  associated  with  relatively 
long  and  straight  fibers. 


THE  STATE  OF  ORIENTATION  OF  THE  FIBER  IS  SPECIFIED 
BY  THE  HERMANS  ORIENTATION  FACTOR f  F* 

FOR  PLANAR  RANDOM  F  =  0 
FOR  PERFECTLY  ALIGNED  F  =  !♦ 

IF  YOU  WISH  DATA  FOR  A  SPECIFIED  ORIENTATION  ENTER  A  "0“  (ZERO) 

IF  YOU  WISH  TO  SCAN  OVER  A  RANGE  OF  ORIENTATIONS  ENTER  A  *1“  (ONE) 

©0 


^  ENTER  SPECIFIC  STATE  OF  ORIENTATION 

©0^ 


HThe  input  data  is  summarized  for  aonvenienoe  and  the 

predicted  properties  displayed  for  the  various  orientations. 

INPUT  DATA 


YOUNG'S  MODULUS 

SHEAR  MODULUS 

POISSON'S 

RESIN 

.5100E+06 

.1960E+06 

.301 

FIBER 

.1050E+08 

.3940E+07 

.333 

FILLER 

♦6930E+07 

.2620E+07 

♦  323 

UOLUME  fractions: 

RESIN  .504 
FIBER  .179 
FILLER  .317 

ASPECT  RATIO  500.0 


CALCULATED  DATA 
ORIENTATION 

LONGITUDINAL  YOUNG'S  MODULUS 
TRANSMERSE  YOUNG'S  MODULUS 
PERPENDICULAR  YOUNG'S  MODULUS 
2y3  SHEAR  MODULUS 

I1.3  SHEAR  MODULUS 
1,2  SHEAR  MODULUS 

1.2  POISSON'S  RATIO 

1.3  POISSON'S  RATIO 

2.3  POISSON'S  RATIO 

STOP 

END  OF  EXECUTION 
CPU  time:  0.77  ELAPSED  TIME:  3:i5.28 
EXIT 


0.00 

.2011E+07 
.2011E+07 
.1615E+07 
.6183E+06 
.6183Ef06 
.7789E+06 
.291 
.287 
♦  287 
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EXAMPLE  2a 


The  purpose  of  this  example  is  to  illustrate  the  use 
of  SMC3  for  a  two  component  system.  In  this  example,  the 
fiber  preperties  will  be  altered  to  reflect  fiber  anisotropy 
(e.g. ,  KEVLAR  49) . 


□  Data  inputs  by  the  user  will  be  indicated  by  the  symbol 


□Comments  are  given  in  italics. 


ENTER  NUMBER  OF  COMPONENTS J 

FOR  RESIN/FIBER  SYSTEM  NUMB  =  2 

FOR  RESIN/FIBER/FILLER  SYSTEM  NUMB  =  3 

TYPICAL  PROPERTIES  FOR  A  TWO  COMPONENT 
POLYESTER/E-GLASS  FIBER  SYSTEM  ARE: 


YOUNG'S  MODULUS  SHEAR  MODULUS  POISSON'S  RATIO 


RESIN  .5100E+06  ♦1960E+06  .301 

FIBER  .1050E+08  .3940E+07  .333 


QThe  stored  values  for  the  fiber  properties  are  to  be  altered. 
The  resin  properties  will  be  maintained. 


IF  THESE  VALUES  ARE  ACCEPTABLE »  ENTER  A  *0"  (ZERO) 

IF  YOU  WISH  TO  USE  SIGNIFICANTLY  DIFFERENT  VALUES  ENTER  A  "1*  (ONE) 

IF  RESIN  PROPERTIES  ARE  ACCEPTABLE  ENTER  A  'O*  (ZERO) 

^  IF  THEY  ARE  TO  BE  ALTERED  ENTER  A  "I*  (ONE) 

Oo 

IF  FIBER  PROPERTIES  ARE  ACCEPTABLE  ENTER  A  “0*  (ZERO) 

^  IF  THEY  ARE  TO  BE  ALTERED  ENTER  A  'I*  (ONE) 

Oi 
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IF  FIBER  IS  ISOTROPIC  ENTER  A  “0"  (ZERO) 

IF  FIBER  IS  ANISOTROPIC  (E*G.  GRAPHITE)  ENTER  A  “1“  (ONE) 
QThe  fiber  is  transversely  isotropic 


Oi 

ENTER  LONGITUDINAL  YOUNG'S  MODULUS 
Ol3«3E6 

ENTER  TRANSVERSE  YOUNG'S  MODULUS 

0:t«33E6 

ENTER  SHEAR  MODULUS)-  G12 

06.88E6 

ENTER  SHEAR  MODULUS »  G23 

Q,6&8£6 

ENTER  POISSON'S  RATIOr  P0S12 

ENTER  POISSON'S  RATIOr  P0S23 

[jThe  new  properties  are  display ed  for  verification  by 
the  user. 


THE  CURRENT  SET  OF  PROPERTIES  FOR  THE  TWO 
COMPONENT  SYSTEM  AREt 


RESIN 

FIBER 

El 

.5100E+06 

.1830E+08 

E2 

.5100E+06 

.1830E+07 

E3 

,5100E+06 

.1830E+07 

G12 

»1960E+06 

.6880E+07 

G13 

♦1960E+06 

.6880E+07 

G23 

.1960E+06 

.688.0Ef06 

P0S12 

♦  301 

♦  300 

P0S13 

.301 

.300 

P0S23 

.301 

.300 

IF  THESE  VALUES  ARE  ACCEPTABLE f  ENTER  A  'O*  (ZERO) 

IF  YOU  WISH  TO  USE  SIGNIFICANTLY  DIFFERENT  VALUES  ENTER  A  "1"  (ONE) 

Oo 
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IF  THE  COMPOSITION  IS  TO  BE  ENTERED  AS  VOLUME  FRACTIONSy 
ENTER  A  “O*  (ZERO) 

IF  THE  COMPOSITION  IS  TO  BE  ENTERED  AS  WEIGHT  FRACTIONS » 
ENTER  A  'I*  (ONE) 

Qo 

VOLUME  FRACTION  RESIN 

VOLUME  FRACTION  FIBER 

0*47 

ENTER  ASPECT  RATIO  OF  THE  FIBER 

O2* 

QThis  low  value  for  the  aspect  ratio  is  associated  with 

pronounced  fiber  curvature  and/or  very  short  fiber  lengths. 


THE  STATE  OF  ORIENTATION  OF  THE  FIBER  IS  SPECIFIED 
BY  THE  HERMANS  ORIENTATION  FACTOR »  F. 

FOR  PLANAR  RANDOM  F  =  0 
FOR  PERFECTLY  ALIGNED  F  =  1* 

IF  YOU  WISH  DATA  FOR  A  SPECIFIED  ORIENTATION  ENTER  A  "0"  (ZERO) 

IF  YOU  WISH  TO  SCAN  OVER  A  RANGE  OF  ORIENTATIONS  ENTER  A  “1“  (ONE) 


'■\Z\The  user  elects  to  scan  over  the  possible  range  of 
orientation. 


ENTER  THE  STARTING  VALUE  FOR  “F" 


ENTER  THE  FINAL  VALUE  FOR  “F* 

ENTER  THE  INCREMENTS  FOR  STEPPING  VALUES  OF  "F" 
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□  input  data  is  summarized  for  oonvenienoe  and  the 
predicted  properties  display ed  for  the  various  orientations . 


INPUT 

DATA 

RESIN 

FIBER 

El 

.5100E+06 

.1830E+08 

E2 

.5100E+06 

.1830E+07 

E3 

.5100E+06 

.1830E+07 

G12 

.1960E+06 

.6880E+07 

G13 

.1960E+06 

.6880E+07 

G23 

.1960E+06 

.6880E+06 

P0S12 

.301 

.300 

P0S13 

.301 

.300 

P0S23 

.301 

.300 

VOLUME 

fractions: 

RESIN  .530 
FIBER  .470 


ASPECT  RATIO  2.0 


CALCULATED  DATA 


ORIENTATION 

LONGITUDINAL  YOUNG'S  MODULUS 
TRANSVERSE  YOUNG'S  MODULUS 
PERPENDICULAR  YOUNG'S  MODULUS 
2f3  SHEAR  MODULUS 

1,3  SHEAR  MODULUS 

1,2  SHEAR  MODULUS 

1.2  POISSON'S  RATIO 

1.3  POISSON'S  RATIO 

2.3  POISSON'S  RATIO 


0.00 

.3065E+07 
.3065E+07 
.1158E+07 
.8269E+06 
.8269E+06 
.1300E+07 
.  178 
.439 
.439 


ORIENTATION 

LONGITUDINAL  YOUNG'S  MODULUS 
TRANSVERSE  YOUNG'S  MODULUS 
PERPENDICULAR  YOUNG'S  MODULUS 

2.3  SHEAR  MODULUS 

1.3  SHEAR  MODULUS 

1,2  SHEAR  MODULUS 

1.2  POISSON'S  RATIO 

1.3  POISSON'S  RATIO 

2.3  POISSON'S  RATIO 


0.25 

.3547E+07 
.2563E+07 
.1155E+07 
.7211Ef06 
.9327E+06 
.1304E+07 
.  215 
.455 
.418 


ORIENTATION 

LONGITUDINAL  YOUNG'S  MODULUS 

TRANSVERSE  YOUNG'S  MODULUS 

PERPENDICULAR  YOUNG'S  MODULUS 

2f3  SHEAR  MODULUS 

1^3  SHEAR  MODULUS 

lf2  SHEAR  MODULUS 

ir2  POISSON'S  RATIO 

1»3  POISSON'S  RATIO 

2»3  POISSON'S  RATIO 


ORIENTATION 

LONGITUDINAL  YOUNG'S  MODULUS 
TRANSVERSE  YOUNG'S  MODULUS 
PERPENDICULAR  YOUNG'S  MODULUS 
2r3  SHEAR  MODULUS 
If 3  SHEAR  MODULUS 
lf2  SHEAR  MODULUS 
if 2  POISSON'S  RATIO 
If 3  POISSON'S  RATIO 
2f3  POISSON'S  RATIO 


ORIENTATION 

LONGITUDINAL  YOUNG'S  MODULUS 
TRANSVERSE  YOUNG'S  MODULUS 
PERPENDICULAR  YOUNG'S  MODULUS 
2f3  SHEAR  MODULUS 
lf3  SHEAR  MODULUS 
If  2  SHEAR  MODULUS 
If 2  POISSON'S  RATIO 
If 3  POISSON'S  RATIO 
2f3  POISSON'S  RATIO 

STOP 

END  OF  EXECUTION 

CPU  time:  0*99  ELAPSED  TIME: 

EXIT 


0.50 

♦4027E+07 
.2065E+07 
.1149E+07 
♦6152E+06 
.1039E+07 
.1300E+07 
♦  266 
.466 
.395 


0.75 

.4505E+07 
.1576E+07 
. 1135E+07 
.5094E+06 
.1144E+07 
. 1285E+07 
.341 
.470 
.373 


1.00 

.4987E+07 
. 1107E+07 
.1107E+07 
.4036E+06 
.1250E+07 
. 1250E+07 
.459 
.  459 
.352 


3:52.92 
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EXAMPLE  2b 


The  purpose  of  this  example  is  to  illustrate  the  use 
of  SMC 3  to  obtain  the  properties  of  a  two  component  ’oni- 
directional  laminate  of  continuous  glass  fibers .  In  this 
example,  stored  properties  will  be  used. 

□Data  input  by  the  user  will  be  indicated  by  the  symbol  O 

□  Comments  are  given  in  italics. 


ENTER  NUMBER  OF  COMPONENTS} 


FOR  RESIN/FIBER  SYSTEM  NUMB  =  2 
FOR  RESIN/FIBER/FILLER  SYSTEM  NUMB  = 


TYPICAL  PROPERTIES  FOR  A  TWO  COMPONENT 
POLYESTER/E-GLASS  FIBER  SYSTEM  AREJ 


YOUNG'S  MODULUS  SHEAR  MODULUS  POISSON'S  RATIO 


RESIN  .5100E+06  .1960E+06  *301 

FIBER  .1050E+08  .3940E+07  *333 


IF  THESE  VALUES  ARE  ACCEPTABLE »  ENTER  A  "0“  (ZERO) 

IF  YOU  WISH  TO  USE  SIGNIFICANTLY  DIFFERENT  VALUES  ENTER  A  "1" 


IF  THE  COMPOSITION  IS  TO  BE  ENTERED  AS  VOLUME  FRACTIONS)- 
ENTER  A  'O*  (ZERO) 

IF  THE  COMPOSITION  IS  TO  BE  ENTERED  AS  WEIGHT  FRACTIONS » 
ENTER  A  (ONE) 

Oo 

VOLUME  FRACTION  RESIN 

VOLUME  FRACTION  FIBER 

0.4 


( ONE ) 
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ENTER  ASPECT  RATIO  OF  THE  FIBER 

Q 1000000. 

\3The  value  of  1000000  fov  the  aspect  patio  is  used  to 
vepvesent  a  continuous  fiber. 


THE  STATE  OF  ORIENTATION  OF  THE  FIBER  IS  SPECIFIED 
BY  THE  HERMANS  ORIENTATION  FACTOR,  F. 

FOR  PLANAR  RANDOM  F  =  0 
FOR  PERFECTLY  ALIGNED  F  =  1. 


OSince  the  properties  of  a  unidirectional  laminate  are 
desired^  the  orientation  factor  is  set  at  unity  to  represent 
aligned  fibers. 


IF  YOU  WISH  DATA  FOR  A  SPECIFIED  ORIENTATION  ENTER  A  "0"  (ZERO) 

IF  YOU  WISH  TO  SCAN  OOER  A  RANGE  OF  ORIENTATIONS  ENTER  A  "1"  (ONE) 

Oo 

^  ENTER  SPECIFIC  STATE  OF  ORIENTATION 

r^The  input  data  is  summarized  for  convenience  and  the 
predicted  properties  displayed.  (Note  that  the  "star 
field"  for  the  aspect  ratio  indicates  a  continuous  fiber 
with  an  infinite  aspect  ratio.) 


INPUT  DATA 


YOUNG'S  MODULUS 

SHEAR  MODULUS 

POISSON'S 

RESIN 

FIBER 

.5100E+06 

♦1050E+08 

.19'60E+06 

.3940E+07 

.301 

.333 

VOLUME 

fractions: 

RESIN  ,600 

FIBER  .400 

ASPECT 

RATIO  ******* 

CALCULATED  DATA 


ORIENTATION 

LONGITUDINAL  YOUNG'S  MODULUo 
TRANSVERSE  YOUNG'S  MODULUS 
PERPENDICULAR  YOUNG'S  MODULUS 
2r3  SHEAR  MODULUS 

1,3  SHEAR  MODULUS 
1,2  SHEAR  MODULUS 

1.2  POISSON'S  RATIO 

1.3  POISSON'S  RATIO 

2.3  POISSON'S  RATIO 


1»00 

,4506E+07 
♦1249E+07 
. 1249E+07 
♦4574E+06 
.5589E+06 
♦5589E+06 
.319 
.319 
♦  365 


STOP 


END  OF  EXECUTION 

CPU  TIMEJ  0.54  ELAPSED  TIME: 


i:i9.02 


46 


CALCULATED  DATA 
ORIENTATION 

LONGITUDINAL  YOUNG'S  MODULUS 
TRANSVERSE  YOUNG'S  MODULUS 
PERPENDICULAR  YOUNG'S  MODULUS 
2y3  SHEAR  MODULUS 

1,3  SHEAR  MODULUS 
1,2  SHEAR  MODULUS 

1.2  POISSON'S  RATIO 

1.3  POISSON'S  RATIO 

2.3  POISSON'S  RATIO 


1.00 

.4S06E+07 
.1249Ei07 
. 1249E+07 
.4574E+06 
.5589E+06 
,5589E+06 
.319 
.319 
.365 


STOP 

END  OF  EXECUTION 

CPU  time;  0.54  ELAPSED  TIME? 


IJ  19.02 
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93100  50  CONTINUE 

93200  CALL  AMATIN  <  CSTAR  »  SL  i- NDUM) 

93300  C - COMPUTE  UPPER  BOUND  FOR  A  PARTICULATE  SYSTEM  <AA-=1) 
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Anticipated  Modifications 

The  program  SMC-3  was  constructed  from  basic  FORTRAN 
statements  in  order  to  facilitate  its  transferability.  The 
program  could  be  made  more  efficient  by  using  statements 
unique  to  specific  computers. 

Currently  SMC-3  is  executed  on  a  DEC-10  system.  For  this 
system,  no  "job"  control  cards  are  required.  Consequently, 
there  are  no  file  or  tape  declaration  statements,  no  calls  to 
specific  compilers,  and  no  memory  or  time  limit  specifications. 
The  following  items  are  summarized  for  the  benefit  of  users 
concerned  with  such  requirements. 

-  Tapes  Declared;  SMC-3  uses  only  input  and  output 

tapes 

-  Compiler:  SMC-3  is  FORTRAN-10  compatible 

-  Memory 

Requirements;  SMC-3  uses  less  than  the  default 

memory  limit  on  the  DEC-10 

-  Time  Limit:  CPU  time  is  usually  less  than  one 

second 

It  was  recognized  that  the  READ  and  WRITE  device  number 
would  vary  with  the  user's  computer  system.  In  order  to 
facilitate  transfer,  SMC-3  has  incorporated  integer  variables 
for  input  and  output  device  numbers .  Both  variables ,  NREAD 
for  input  device,  NWRITE  for  output  device,  are  assigned  values 
using  a  DATA  statement  (line  12100  of  the  program) .  The 
typical  input/output  commands  in  SMC-3  are  of  the  form 

DATA  NREAD, NWRITE/ 5, 5/ 

m 

READ (NREAD, 10) A 


WRITE (NWRITE, 100) A 
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To  assign  the  correct  device  numbers  to  the  READ  and  WRITE 
statements,  it  is  only  necessary  to  change  the  DATA  statement 
(line  12100).  The  variables  NREAD  and  NWRITE  are  transferred 
to  the  required  subroutines  through  COMMON/Bl/.  For  ease  of 
data  entry,  SMC-3  was  written  using  free  formatted  READ 
statements.  The  free  format  READ  symbol  for  the  DEC-10  is 
the  star,  These  READ  statements  are  of  the  form: 

READ (NREAD, *)B 

If  the  computing  system  on  which  SMC-3  is  being  implemented 
has  free  format  READ  capabilities,  the  correct  symbol  will  need 
to  be  used  in  place  of  the  star.  For  systems  which  do  not  have 
the  free  format  READ  options,  it  will  be  necessary  to  format 
each  of  the  existing  free  format  READ  commands.  Young's  and 
Shear  moduli  could  be  read  using  the  exponential  field  format — an 
E10.5  field  would  suffice.  All  other  formats  could  be  replaced 
by  a  floating  point  field.  For  example,  a  F20.10  would  suffice. 
If  this  is  done,  all  data  should  be  entered  including  decimal 
points .  Statement  numbers  for  the  added  FORMAT  statements 
should  be  500  or  longer  to  prevent  any  duplication  of  existing 
statement  numbers . 

To  replace  the  free  formatted  READ  statements,  the  following 
changes  will  be  necessary: 

1)  change  free  format  symbol,  to  FORMAT  statement 

numbers 

2)  insert  corresponding  FORMAT  statement 

For  example,  the  READ  statement  which  reads  the  volume  fractions 
(program  line  30300)  could  be  changed  to: 

READ (NREAD, 500) V2(K) 

500  FORMAT (F20. 10) 

The  READ  command  for  Young's  modulus  (program  line  46500)  could 
be  rewritten  as: 

READ (NREAD, 510)  El  (I) 

510  FORMAT  (ElO.5) 


A  listing  of  each  line  where  the  free  format  has  been 
utilized  is  summarized  below.  Lines  30300  through  36800  are 
all  contained  in  SUBROUTINE  INPUT.  All  these  lines  can  be 
replaced  with  a  single  floating  point  FORMAT  statement.  The 
remaining  free  formatted  READ  statements  occur  in  SUBROUTINE 
ALTER.  Both  exponential  field  and  floating  point  field  FORMAT 
statements  will  need  to  be  used  in  this  siibroutine.  POS12(2), 
POS23(2)  and  POS  should  be  read  using  the  floating  point  field 
while  the  remaining  variables  can  be  read  using  the  suggested 
exponential  field. 


(From  SUBROUTINE  INPUT) 

30300 

READ(NREAD»*) 

02<K> 

31800 

REAIUNREADf*) 

DX  ( K ) 

32100 

"  REAIKNREADr*) 

DEN<K) 

34300 

REAIKNREADf*) 

AA2 

35600 

READ(NREAD»*) 

FSTART 

36200 

REAIKNREADf*) 

FSTART 

36500  - 

REAEKNREADf*) 

FSTOP 

36800 

REAEKNREADr*) 

FADEi 

(From  SUBROUTINE  ALTER) 

44200 

REAEKNREAEif*) 

El(2) 

44500 

REAEKNREAEi!-*) 

E2(2) 

44900 

READ(NREAD»*) 

G12(2) 

45300 

REAEKNREADi-*) 

G23(2) 

45600 

REAEi<NREAEi»*) 

P0S12<2) 

46000 

REAEKNREAEi)-*) 

P0S23(2) 

46500 

REAEKNREAEi)-*) 

E1<I> 

47000 

REAEKNREAEif*) 

G12<I) 

47500 

REAEKNREAEif*) 

POS 

* 


Another  change  which  may  be  necessary  deals  with  contin¬ 
uation  cards.  In  two  FORMAT  statements,  the  number  of 
continuation  cards  has  exceeded  four.  Seven  were  used  starting 
at  line  20000  while  six  were  used  starting  at  line  34600.  For 
computers/compilers  which  are  limited  to  fewer  continuation 
cards,  the  output  at  these  two  locations  will  need  to  be 
rewritten  using  two  WRITE  statements. 


SPECIALIZED  TI-59  ROUTINES 


Introduction 


This  section  describes  the  operation  of  a  TI-59  calculator/ 
PC-100  printer  programmed  to  predict  properties  for  two-component 
(fiber/resin)  and  three  component  (fiber/filler /resin)  sheet 
molding  materials. 

The  program  is  segmented  on  four  magnetic  cards: 

Card  I  —  Reads  input  and  generates  reference  phase 
Card  II  —  Generates  CSTAR 
Card  III  —  Planar  averaging 

Card  IV  —  Generates  Engineering  constants  and 
controls  output. 

The  procedures  for  reading  magnetic  cards  are  reviewed  at  the 
end  of  this  section. 

A  PC-100  printer  is  required  for  the  operation  of  the 
program. 

The  current  version  of  the  program  is  restricted  to  iso¬ 
tropic  fibers  with  aspect  ratios  in  excess  of  150.  Supplementary- 
cards  will  be  provided  at  a  future  date  to  deal  with  low  aspect 
ratio  fibers  and  platelet  reinforcing  agents . 

The  following  input  data  is  required: 


Ej  =  Young's  modulus  of  component  "j" 


Vj  =  Poisson's  ratio  of  component  "j' 


Gj  =  Shear  modulus  of  component 


V.  =  voliame  fraction  of  component 


j 


f  =  orientation  parameter,  0  £  f  £  1 

(f=0  is  random,  f=l  is  perfectly  aligned) 

The  operating  procedures  for  two-phase  and  three-phase  systems 
are  summarized  in  the  following  section.  Sample  calculations 
are  given  in  the  subsequent  section. 


Preprogrammed  magnetic  cards  will  be  provided.  However, 
difficulty  has  been  encountered  in  reading  magnetic  cards 
programmed  on  other  machines.  Consequently,  a  program  listing 
is  provided  so  that  cards  can  be  generated. 


□  Card  II 

□  Read  sides  1  and  2  of  card  II 

□  Press  A 

Printer  will  respond  with  a  print 
and  advance  when  finished 

□  Planar  averaging 

□  Read  sides  1  and  2  of  card  III 
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Press 

□  Enter  orientation  parameter  f  A 

Planar  averaging  B 

□  Output 

□  Read  sides  1  only  of  Card  IV 

□  Generate  output  A 


^SUMMARY  OF  RESULTS* 


I>To  vary  f,  the  orientation  parameter: 
Run  the  program  as  before,  with  the 
first  value  of  f. 


> After  obtaining  the  necessary  output, 
read  side  1  (only)  of  card  III 

>Enter  orientation  parameter,  f,  and 
continue  with  planar  averaging 
operation 
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OPERATING  PROCEDURES  FOR  THREE-PHASE  SYSTEMS 


□  Enter  physical  properties  of  resin  and  filler 
□  Read  sides  1  and  2  of  card  I 


□  Enter  resin  properties 


□  Enter  filler  properties 


°=Fill 


°°Fill 


Press 


Printer 


A 

^R 

S 

B 

V 

R 

S 

C 

S 

2nd 

A' 

^Fill 

Sill 

2nd 

B' 

^Fill 

^Fill 

2nd 

C 

Sill 

Sill 

□  Enter  modified  volume 
fraction  filler 

,  ^  '^Flll 

Fill 

□  Generate  surrogate  matrix 
properties 


^Fill 


^Fill 


□  Enter  physical  properties  of  surrogate 
matrix  and  fiber 

□  Enter  surrogate  matrix  properties 
□  E 


C 
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REVIEW  OF  CARD  READING  PROCEDURES 
I .  Reading  Cards 

All  program  cards  need  both  sides  read  except  the  output 
program.  When  an  instruction  tells  you  to  read  both  sides  of  a 
card: 

A.  Press  CLR 

B.  Slide  card  through  lower  slot  in  calculator,  printed 
side  up: 


The  card- reader  motor  will  start.  Push  the  card  in  until 
the  gears  catch  it  and  pull  it  through.  Pull  the  card 
out  from  the  other  side.  A  1  should  be  in  the  display. 

C .  Press  CLR 

D.  Turn  card  around  and  slide  through  again  (printed  side 
still  up. 


E.  If  display  flashes,  press  CLR  and  try  again. 

F.  If  motor  continues  to  run  and  display  is  blank  after 
reading  a  card,  press  R/S .  There  wasn't  anything  on 
the  side  you  read.  Make  sure  the  card  was  supposed  to 
have  both  sides  read. 
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EXAMPLES 
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EXAMPLES  RUN  WITH  A  TWO-PHASE  COMPOSITE 

SMC- 6 5  Properties  (in  psi) 


E  V  G  Vol .  Fraction 


Resin  5.1  x  10^  .301 

Fiber  1.05  x  10^  .333 

1.96  X  10^ 
3.94  X  10® 

.6 

.4 

□  Read  sides  1  and  2  of  card  I 

Press 

Di^ 

Printer 

Enter  resin  properties 

□  5.1  EE  5 

A 

5.1 

05 

□  .301 

B 

3.01 

-01 

3.01-01  ’  '  ' 

□  1.96  EE  5 

c 

1.96 

05 

1 i 0^ 

Enter  fiber  properties 

□1.05  EE7 

2nd  A' 

1.05 

07 

1.05  07  :  . 

□  .333 

2nd  B' 

3.33 

-01 

3=  33  -0 1 

□3.94  EE  6 

2nd  C 

3.94 

06 

3=  94  06 

'  ... 

Volume  fraction  fiber 

□  .4 

D 

4. 

-01 

■  4. 

□  Generate  reference  phase 

. 

properties 

E 

0 

s  ^  *».* 

2. 821-01 

□  Read  sides  1  and  2  of  card  II 

5.  225  05 

□  Run  program  on  this  card 

A 

0. 

s 

□  Read  sides  1  and  2  of  card  III 
Enter  f 

□  0.  AO 

□  Planar  averaging  B  0  Printer  advance 

when  done 
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□  Read  side  1  of  card  IV 

□  Generate  output 

♦SUMMARY  OF  RESULTS* 


□  Read  side  1  of  card  III 

Enter  new  f  value 
□  .25 

□Planar  averaging 

□  Read  side  1  of  card  IV 

□Generate  output 

♦SUMMARY  OF  RESULTS* 
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SAMPLE  RUN  WITH  A  THREE-PHASE  COMPOSITE 


SMC-25  Properties  (in  psi) 


E 

V 

G 

Vol.  Fraction 

Resin 

5.1  X  10^ 

.301 

1.96  X  10^ 

.504 

Fiber 

1.05  X  10^ 

.333 

3.94  X  10^ 

.179 

Filler 

3  X  10^ 

.402 

1.07  X  10® 

.317 

□Read  sides  1  and  2  of  card  I 

Enter  resin  properties 

□  5.1  EES 

□  .301 

□  1.96  EES 

Enter  filler  properties 

□  3  EE6 

□  .042 

□  1.07  EE  6 

Enter  modified  volume  fra 
filler; 

V  .  _ 

''Fill  ''r 


°''f111  °  ^  .504  °  -3861  D  3.861-01 

□  Generate  surrogate  matrix 
properties 

Enter  surrogate  matrix 
properties 

□9.582  EES 

□  .3117 

□  3.652  EES 

Enter  fiber  properties 

□1.05  EE  7 

□  .333 

□  3.94  EE  6 


A  9.582  05 
B  3.117-01 
C  3.652  05 


2nd  A'  1.05  07 

2nd  B'  3.33  -01 
2nd  C  3.94  06 


Press 


Printer 


Volume  fraction  fiber 
□  .179 


D  1.79  -01 


□  Generate  reference  phase 
properties 


□  Read  sides  1  and  2  of  card  II 

□  Run  program  on  this  card  A  0.  0;.  v 

□  Read  sides  1  and  2  of  card  III 

Enter  f 

□  0  AO  U . 

□  Planar  averaging  BO  1  advance 


□  Read  side  1  of  card  IV 
□  Generate  output 


♦SUMMARY  OF  RESULTS* 
(f=0) 


□  Read  side  1  of  card  III 

Enter  new  f  value 
□  .25 

□  Planar  averaging 

□  Read  side  1  of  card  IV 

□  Generate  output 

♦SUMMARY  OF  RESULTS* 
(f=0.25) 


TI-59 

PROGRAM  LISTINGS 

The  key  strokes  associated  with  cards  I-IV  are 
listed  on  the  following  pages . 

It  is  recommended  that  upon  entering  the  routines, 
duplicate  magnetic  cards  should  be  recorded.  The  spare 
card  serves  as  a  replacement  for  a  damaged  card.  If  a 
card  is  damaged,  a  second  duplicate  should  be  recorded 
from  the  spare  card. 
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CARD  I 
-continued- 


rv 

03  3 

“j  ••  *1  -“i 

O  w 

71  SBR 

-Jl*  •  .1 

'i-2 

STD 

274 

44  SUM 

*  "l  ‘"^1 

^*2  O 

6  3  H  0  P 

-“*1  -**1 
vZ'  (  iC 

10 

10 

V  /■  ~i 

00  00 

324 

61  GTO 

A  -*“» 

'i'  O 

RCL 

2  ^"6 

76  LBL 

325 

78  It 

374 

07 

07 

iC  j  r‘ 

3  4  r  A 

326 

76  LBL 

375 

J  .—I 

STD 

•-I  i“i 

i'  o 

03  3 

1 

68  HDP 

3  7  6 

1 1 

1 1 

279 

■74  SM* 

328 

43  RCL 

7 1 

SBR 

280 

09  09 

3  2  9 

4  b  4 

i**i 

iL.  tL. 

IHV 

281 

97  IiS2 

330 

64  Plit 

379 

1  O 

RCt 

*1*  C* 

Cm  im. 

09  09 

331 

06  06 

380 

00 

00 

•'**1 

.■*  Cl*  O 

34  TK 

o  o  ^ 

6  4  P  D 

C:  1 

1^1  ^ 

74 

SMt 

284 

71  SBR 

-1  o 

07  07 

382 

06 

06 

285 

52  EE 

334 

6  4  P  D 

*Z^  i^i  o 

('  b 

R  C  ^ 

23  6 

01  1 

335 

08  08 

384 

01 

01 

tC*  J  1 

94  +/- 

336 

43  RCL 

385 

74 

SM# 

L^!  c* 

44  SUM 

3  'I*  f  * 

43  43 

3  8  6 

07 

07 

2  Q  9 

42  42 

“1  •  •»  o 

O  *—* 

55  t 

•T‘  O  ”7 

71 

SBR 

280 

43  RCL 

!  "1 9 

4  3  R  L-  L 

■Z'  O  b 

t  mm 

INV 

291 

42  42 

340 

42  42 

389 

92 

RTN 

292 

35  1/K 

341 

95  ^ 

390 

i‘  b 

LBL 

293 

85  + 

3  *-7  2 

**?  A  t~’  M  ••• 

|-  ^ 

39 1 

i"  O 

Z-i- 

294 

01  1 

•~i  A  •“» 

0  6  0  6 

3  2 

01 

1 

295 

95  = 

344 

43  RCL 

393 

75 

- 

*P  Q 1^, 

42  STu 

A  CT 

4^1  4^x 

44  44 

3  9  4 

43 

RCL 

297 

4  6  4  6 

3  4  6 

55  t 

395 

•H-  2 

'A-  A 

'"i  J-j 

7’’  1  3  R 

3  i-  r' 

43  RCL 

'396 

95 

- 

299 

68  HOP 

3  *^1’  o 

42  42 

397 

42 

STO 

300 

01  1 

3  4  9 

95  = 

' Q  !-! 

46 

46 

3  L»  1 

44  SUM 

350 

74  SM* 

399 

cr  nr. 

“T 

302 

42  42 

351 

07  07 

400 

02 

2 

303 

4  3  R  L-  L 

352 

43  RCL 

401 

65 

3  ^  j  *4 

4  b  •t  b 

3  5  3 

45  45 

402 

43 

RCL 

305 

35  1/K 

354 

55  t 

A  1 

“i-UO 

42 

4  2 

306 

H  -f  c  T  n 
‘t  L_  1  la 

355 

4  3  RCL 

.-i  t”:  A 
i  i 

95 

= 

307 

46  4 €• 

356 

42  42 

OR  4. 

4U5 

42 

A 

ri 

•I«  1  U 

3  w 

309 

42  sfn 

3  5  7 
358 

O  •  J  ^ 

73  RC^ 

*T  »— •  O 

4  U  7 

03 

t" 

310 

09  09 

359 

08  '  08 

408 

A  'Z* 

-r 

■Zt  i  U 

3 1 1 

03  3 

360 

95  = 

409 

09 

09 

‘  1  2 

94  +..•••- 

36 1 

35  1--"X 

410 

0 1 

1 

313 

44  SUM 

3  E*  2 

85  H- 

4 1  1 

4  n' 

•z- 1  i  hi 

314 

00  00 

363 

-?  r’t  »“•  •. 

r  o  r.L.T- 

412 

08 

0  8 

31 5 

76  LBL 

364 

02  02 

413 

76 

LBL 

3 1 6 

89  -n’ 

365 

95  = 

414 

70 

RHil 

• 

74  SM* 

3  6  6 

35  1/K 

415 

4  3 

RCL 

O  O 

09  09 

367 

7  9  T  i 

416 

46 

4  b 

319 

97  iiSZ 

368 

08  08 

4 1  7 

65 

320 

09  09 

369 

43  RCL 

4 1 3 

b 

RC4 

3  2 1 

89  iT 

370 

n  6  u  fct 

419 

U  E‘ 

f\f, 

\/ 

{/ 

96 


CARD  I 
-continued- 


NOTE:  Intermediate  results  for  CZERO  are  stored  in  the 

following  regist(Ers 


Reg  36 : 

^11  ^22  ^33 

37; 

^12  ^  ^13  ^  ^23 

^44  “  ^55  ^66 


38 
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CARD  II 
-continued- 


135 

35  1/K 

184 

42  STD 

2  -L*  3 

85 

4 

136 

04  4 

135 

12  12 

234 

43 

RCL 

137 

U  6  6 

186 

43  RCL 

235 

26 

•^1  ^ 

138 

42  STQ 

187 

43  43 

2  •-*  b 

95 

= 

139 

05  05 

188 

85  + 

2*:>  ( 

35 

140 

06  6 

189 

43  RCL 

•“«  o 
.z« 

85 

4 

141 

'42  STO 

190 

•“i  ”“i  •"«  ’Ij 

end  side  1-239 

43 

RCL 

1  4Ic! 

00  00 

191 

95  = 

240 

•‘^1  o 

'“i  o 
1^1 

143 

01  1 

192 

42  STD 

24  1 

95 

= 

1 44 

0  6  b 

193 

13  13 

b  4  b 

42 

STD 

cr 

i  ‘-r.J 

42  STO 

194 

43  RCL 

2  *■'+  3 

4  6 

A 

••t  o 

t 

07  07 

195 

t::|.  <4  4 

244 

25 

CLR 

1 47 

76  LBL 

196 

85  + 

245 

99 

PRT 

1 43 

49  PRD 

197 

43  RCL 

246 

98 

flDV 

i  4  y 

73  RC* 

198 

24  24 

24  i*' 

9 1 

R  ♦  3 

150 

05  05 

199 

95  = 

7  b 

LBL 

151 

65  X 

200 

42  STD 

2'-t  9 

■"1  .-i 

il  •i’ 

Cl 

152 

43  RCL 

201 

14  14 

250 

03 

1  ’J  O 

40  40 

202 

f  1  5u'K 

•-J  nr  T 
il-J  X 

42 

i  1  1 

i  LJ 

154 

85 

203 

35  1/K 

2b  2 

00 

00 

155 

73  RC^ 

2U4 

43  RCL 

J  Sb* 

faZ*  •— *  O 

i  •— * 

Lb  L 

156 

06  06 

205 

3  6  3  6 

254 

53 

( 

1 57 

6  5 

2  0  6 

44  sun 

255 

4  3 

R  L-  L 

158 

43  RCL 

207 

41  41 

256 

04 

04 

159 

50  50 

2  o 

44  sun 

2b  7 

r^b 

- 

160 

95  = 

209 

43  43 

258 

01 

•i 

X 

1  i:  1 

!*  2  5  1 

2 1 0 

.■i  •'“i  r"t  ! 

o  r-.  L-  L- 

2  b  9 

95 

= 

162 

07  07 

211 

O  i'  1'  1 

260 

*7  O 

l  w 

ST-^ 

1 63 

72  ST* 

•-1  i  O 
tl.  4 

44  sun 

ill  O  X 

00 

00 

164 

0 5  0  5 

•-*1  o 

il.  X  •-.* 

42  42 

2  o  2 

4  ki 

r-  -r  n 
O  1  U 

165 

01  1 

-i  A 
il  4  •i- 

44  sun 

2  o 

04 

s“i  .•! 
tj  «-r 

166 

94 

O  1  cr 

xX,  4  *..> 

44  44 

•"t  . "  A 

rZ  O  ^ 

97 

4 

A  ‘“i  r"l  1 

“J/’ 

r* 

f  1  f\ 

1  !* 

oUrl 

il  1  o 

H*  -I*  r.  »-•  L 

iZ  C‘  •«.» 

1^:  s^t 

168 

05  05 

217 

45  45 

2  6  6 

53 

(. 

1  »r‘9 

44  sun 

35  1/K 

ilL  C'  i' 

R*  1 !  4 

170 

06  06 

2 1 9 

85  + 

2  63 

01 

01 

44  sun 

220 

4  3  R  C  L 

2  9 

b  5 

- 

1 

X  > 

07  07 

hi  ii~>  1 

25  25 

270 

4  -1* 

RCL 

1  ^3 

97  DSZ 

O  *1*  o 

95  = 

0  7  i 

36 

36 

1  Td 

i  < 

00  00 

•;:*  c* 

w  k— 

35  1/K 

b  *7 

95 

= 

175 

49  PRB 

224 

85  H- 

4  '  3 

42 

STO 

176 

71  SBR 

225 

43  RCL 

274 

10 

10 

177 

35  1/K 

226 

O  *1*  ‘Ij  jlj 

.-.-pc: 

y  'ij 

RC^ 

178 

43  RCL 

.-I  ■:i  “7 
i4.  i4.  i 

95  = 

276 

02 

02 

179 

41  41 

*1‘  *  -*  O 

^  iX. 

42  STD 

•”“t  ■*?**? 

li.  r 

75 

- 

180 

,•1  I*"*  "T  n 

•T;C.  •!«  1  U 

229 

4b  45 

278 

A  ♦“I 

•i*  C‘ 

RCL 

i  C«  i 

X  i.J  X 

11  11 

230 

43  RCL 

bbQ 

O  "7 

— ;i 

182 

43  RCL 

O  i 

46  46 

230 

95 

= 

•t  S^X 

sr 

“i  :“j 

.■5  “I 

C’  T  n 

•Til  ■-i-kl 

*1.  O  2 

O  ‘J  X  '■  « "i 

X 

_«  1  kJ 
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CARD  II 
-continued- 


.-1  r: 

l"j 

i  1 

1 

331 

43 

RCL 

380 

43 

RCL 

O  C’  O 
iZ.  ■— * 

71 

SBR 

.  C*  ilX 

13 

1 3 

381 

1  1 

1 1 

d  4 

Cmm  %Lm 

INV 

333 

75 

- 

3  0  2 

75 

- 

2S5 

-?  --J 

l‘  0 

RC* 

334 

43 

RCL 

3  3  3 

43 

RCL 

286 

03 

03 

335 

Cm 

25 

384 

12 

•i  •“1 

1  cZ 

-1  o  -7 
iZw  1 

75 

- 

336 

95 

= 

385 

‘”1  c* 
<Z» 

A- 

.-i  I'l  0 

^  0  * 

id 

RCL 

.-I -y 

J 

•^sc: 

■— * 

i 

i  •  II 

3 b 

42 

STD 

239 

3  y 

’“i  Cl 

0  *— * 

•-I  •-!  i-i 

0  ’Z*  0 

72 

ST* 

*7 

-I*  c»  r 

08 

:”t  i“i 

U  0 

290 

95 

= 

339 

05 

05 

•^1 

w  w 

65 

291 

35 

1/K 

340 

43 

RCL 

389 

02 

c> 

292 

42 

STD 

341 

1  0 

i.  •-» 

4  ■-1 

X  -I* 

390 

95 

= 

«-i  0 
.7  0 

13 

13 

342 

75 

- 

391 

42 

sin 

294 

06 

6 

343 

43 

RCL 

Cm 

09 

09 

295 

42 

STO 

3  4  4 

26 

26 

•  ^7 

35 

1  •**’  K 

296 

00 

00 

345 

95 

= 

! "!  *4 

65 

>- 

297 

7  6 

LBL 

346 

35 

1/X 

395 

43 

RCL 

298 

54 

o4  { 

\  lU 

ST* 

396 

07 

07 

--  «* 

43 

RCL 

31 4  8 

06 

6 

397 

95 

= 

300 

07 

07 

349 

92 

RTH 

398 

"7  '*!’ 

ST4 

301 

75 

- 

350 

7  6 

LBL 

'  -j  ‘-7  M 

01 

01 

302 

01 

■j 

351 

^  J**"! 

f  »«* 

1  + 

400 

it  3 

RCL 

303 

95 

r: 

352 

*7 

1*  0 

RC* 

401 

4  c« 

1  w 

1  2 

304 

-7 

1  L- 

ST* 

353 

01 

01 

*i*  U  iC 

94 

-H/  - 

305 

00 

00 

354 

42 

STO 

403 

55 

-r 

306 

42 

STO 

355 

1 1 

4  4 

i  1 

404 

.4 

RCL 

307 

07 

07 

3!5^6 

r’  3 

RC* 

405 

09 

09 

308 

97 

nsz 

357 

02 

02 

406 

95 

= 

309 

00 

00 

0  0 

42 

STO 

407 

r*  ill 

3T^ 

310 

54 

*) 

359 

12 

12 

408 

U2 

uc. 

3  :  1 

43 

RCL 

360 

-7  .~j 

f  0 

RC* 

409 

^  0 

K  U  L 

0  X  Im 

10 

10 

361 

03 

03 

410 

1 1 

1  1 

,-}  X  •"< 

(  tZ, 

yl  * 

362 

42 

STO 

411 

65 

314 

01 

01 

363 

13 

13 

412 

4  3 

RCL 

3 1 5 

75 

- 

364 

3 

RC* 

413 

i  3 

I  3 

6!  1 6 

43 

RCL 

365 

04 

04 

414 

—*tT 

i  •  J 

- 

3  i  r 

0  " 

Urn  * 

.♦  .^1 

3  6  6 

.4  •“! 
•i-ul 

STO 

415 

43 

RCL 

318 

95 

= 

3  b  7 

14 

14 

416 

08 

08 

3 1 9 

1*’  t:! 

ST* 

3  b  C' 

92 

RTH 

4  i  r 

95 

= 

320 

03 

03 

369 

7  6 

LBL 

413 

55 

4 

321 

43 

RCL 

0  ‘  J 

35 

1/X 

4  1  9 

53 

( 

•~S  ♦“*«  0 
w 

11 

11 

371 

43 

RCL 

420 

53 

( 

0  ^  0 

-7  .-1 

I*  iZ. 

ST* 

•*":  *"’1 

0  (  C. 

13 

13 

A  •’*■«  '5 

*i'c-  i 

4  3 

RCL 

324 

02 

02 

*■?  •“« 

0  i  •-> 

85 

M 

4  id  2 

1 3 

i  *1' 

3  2  5 

75 

- 

374 

43 

RCL 

.4  ‘'j  c» 

- 

.326 

43 

RCL 

375 

14 

14 

424 

43 

RCL 

327 

24 

•“t  .4 

Ll'i* 

376 

95 

= 

425 

•I  ^ 

1  id 

*^4  l"j 

iZ. 

95 

=: 

0  t  ( 

42 

STD 

A  C 

54 

329 

*7  0 

CT:u 

0  \  * 

•“1  *r«  ri 

C>  ;*  *_* 

07 

07 

427 

65 

330 

04 

04 

379 

65 

4::10 

43 

RCL 
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CARD  II 
-continued- 


429 

09  09 

430 

54  ) 

A 

^  .'l  1 

42  STQ 

432 

09  09 

433 

95  = 

434 

T*  ■“’i  1*^  “T  V 

435 

03  03 

436 

43  RCL 

I* 

08  08 

4  o  o 

75  - 

439 

43  RCL 

440 

n  11 

•H"  1 

65  X 

442 

43  RCL 

A  A  •~i 

O 

4  .i  -i  .4 

i  ‘i-  i 

444 

95  = 

445 

55 

4  4 

43  RCL 

44? 

09  09 

448 

95  = 

449 

1*’  2  b  1 

450 

U  4  U  4 

451 

92  RTN 

NOTE:  Intermediate  results  are  stored  in  the  following 

registers. 


EZERO 

CSTAR 

Reg  2 1 ; 

^11 

=^22=“ 

Reg  41; 

^11 

22; 

®12 

=  ®13  '  “ 

42: 

^12 

=  ^13 

23: 

^33 

43; 

^22 

=  S3 

24 

^23 

44: 

^23 

25 

^44 

45: 

^44 

26 

^55 

=  ®66 

46; 

^55 

=  S6 

101 


CARD  III 

Input  Orientation  Factor  Cf)  and  Planar  Averaging 


Labels 

Used 

042 

85  t 

091 

95 

= 

001 

11 

fl 

043 

43  RCL 

092 

42 

STO 

\J  i.  ■H* 

12 

B 

044 

17  17 

093 

09 

09 

£66 

INV 

045 

95  = 

094 

53 

i;| 

046 

42  STO 

095 

A  •-! 

’-r-o 

RCL 

U  r 

19  19 

096 

42 

42 
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